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ABSTRACT 

We present the small-scale (0.2/i~^ to 7h~^ Mpc) cross-correlations between 32,000 luminous early- 
type galaxies and a reference sample of 16 million normal galaxies from the Sloan Digital Sky Survey. 
Our method allows us to construct the spherically averaged, real-space cross-correlation function between 
the spectroscopic LRG sample and galaxies from the SDSS imaging. We report the cross-correlation as a 
function of scale, luminosity, and redshift. We find very strong luminosity dependences in the clustering 
amplitudes, up to a factor of 4 over a factor of 4 in luminosity, and measure this dependence with 
high signal-to-noise ratio. The luminosity dependence of bias is found to depend on scale, with more 
variation on smaller scales. The clustering as a function of scale is not a power law, but instead has a dip 
at lh~^ Mpc and an excess on small scales. The fraction of red galaxies within the L* sample surrounding 
LRGs is a strong function of scale, as expected. However, the fraction of red galaxies evolves in redshift 
similarly on small and large scales, suggesting that cluster and field populations are changing in the 
same manner. The results highlight the advantage on small scales of using cross-correlation methods as 
a means of avoiding shot noise in samples of rare galaxies. 



Subject headings: cosmology: observations — galaxies: clusters: general 
lenticular, cD — galaxies: evolution 



galaxies: elliptical and 



1. INTRODUCTION 

Galaxy clustering allows us to study the relation of galax- 
ies to dark matter through the biased clustering of dark 
matter halos and to probe the possibility of environment- 
dependent processes in the evolution of galaxies. With 
today's large galaxy surveys, we can measure this clus- 
tering at very high signal-to-noise ratio, investigating the 
detailed trends of clustering with luminosity, color, mor- 
phology, and redshift (Hubble 1936; Zwicky et al. 1968; 
Davis & Geller 1976; Dressier 1980; Postman & Geller 
1984; Hamihon 1988; White, TuUy, & Davis 1988; Park 
et al. 1994; Loveday et al. 1995; Guzzo et al. 1997; Benoist 
et al. 1998; Willmer, da Costa & Pellegrini 1998; Brown, 
Webster & Boyle 2000; Carlberg et al. 2001; Norberg et al. 
2001; Zehavi et al. 2002; Norberg et al. 2002; Budavari et 
al. 2003; Madgwick et al. 2003; Hogg et al. 2003a; Zehavi 
et al. 2004b). 

Massive galaxies are particularly interesting to study via 
clustering analyses because they tend to reside in mas- 
sive dark matter halos (Sandage 1972; Hoessel et al. 1980; 
Schneider et al. 1983; Postman & Lauer 1995). These halos 
themselves have strong trends of clustering amplitude ver- 
sus mass (Kaiser 1984; Bardeen et al. 1986; Mo & White 
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1996; Sheth, Mo, & Tormen 2001), implying that one can 
find significant variations in clustering whenever the mean 
halo mass of the subsample galaxies is changed. Moreover, 
because massive galaxies are roughly passively evolving at 
low redshifts (McCarthy et al. 2001; Fontana et al. 2004; 
Glazebrook et al. 2004, but see Drory et al. 2004), we 
should be able to interpret the redshift evolution of clus- 
tering. A fixed set of galaxies cannot be in high-mass halos 
at one redshift and low-mass halos at the next, whereas a 
set defined by a more transitory identifier, such as strong 
star formation, need not be found in a consistent location 
from one time to another. 

Here, we use a sample of very luminous early-type galax- 
ies from the Sloan Digital Sky Survey (SDSS; York et 
al. 2000) to measure the small-scale clustering (0.2ft,^^ 
to 7ft.^^Mpc) of the most massive galaxies. We do this 
by use of angular cross-correlations between 32,000 lumi- 
nous red galaxies (LRG) with spectroscopic redshifts and 
16 million fainter galaxies from the SDSS imaging survey. 
The method of Eisenstein (2003) constructs spherically 
symmetric measurements of the real-space spatial cross- 
correlation function in a flexible manner that supports in- 
vestigations of many different potential trends. We have 
previously applied the method to the study of low-redshift 
galaxies in the SDSS (Hogg et al. 2003a; Blanton et al. 
2003c). We focus here on the luminosity, scale, and color 
dependence of the clustering of ^ L* galaxies surrounding 
the LRGs. 

The application to LRGs shows the considerable advan- 
tages of the cross-correlation of spectroscopic and imaging 
data sets for the study of small-scale clustering of rare pop- 
ulations. By avoiding auto-correlation, we remove a large 
amount of shot noise from the computation at only a minor 
cost in signal. Moreover the flexibility of the method rec- 
ommends it to topics involving trends in multi-dimensional 
parameter spaces. We expect that the method will find 
applications in higher redshift work as well. 
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As tracers of high-density regions, massive early- type 
galaxies also offer a way to study the environmental de- 
pendences of nearby fainter galaxies via cross-correlation 
techniques. We do this by using the colors of the galaxies 
in the imaging sample to measure the fraction of red galax- 
ies as a function of redshift and distance from the LRG. 
Of course, there is a strong gradient in red fraction as a 
function of distance (Hubble & Humason 1931; Abell 1965; 
Oemler 1974; Melnick & Sargent 1977; Dressier 1980; Post- 
man & GcUcr 1984). The time cvohition of the red fraction 
as a function of environment is one of the more studied 
quantities in extragalactic astrophysics, in clusters as the 
Butcher-Oemlcr effect (e.g.. Butcher & Oemler 1978, 1984; 
Dressier & Gunn 1983; Lavcry & Henry 1986) and in the 
field as the faint blue galaxy problem (e.g.. Broadhurst ct 
al. 1988; Tyson 1988; Jones et al. 1991; Lilly et al. 1991; 
Metcalfe et al. 1991) and the evolution of the star forma- 
tion density (e.g., Lilly ct al. 1996; Madau et al. 1996). 
While all agree that these reddening trends are a window 
into galaxy evolution, the questions of whether this evolu- 
tion is environmentally dependent and what processes un- 
derlie it remain controversial (e.g., Hashimoto et al. 1998; 
Balogh et al. 1999; Margoniner et al. 2001; de Propris et 
al. 2004). 

The outline of this paper is as follows. §2 reviews the 
methodology employed, with more detail provided in Ap- 
pendices A and B. §3 describes the samples of spectro- 
scopic and imaging galaxies we use. We present the re- 
sults in §4 and conclude in §5. We use a Qm = 0.3 flat 
cosmology throughout and quote magnitudes with h = 1. 
All magnitudes have been corrected for reddening using 
the Schlegcl et al. (1998) map. 

2. CROSS-CORRELATION METHODOLOGY 

We will use the angular cross-correlation method de- 
scribed in Eisenstein (2003). In this method, one cross- 
correlates a sample of known redshift with a sample of 
unknown redshift. The first sample will be spectroscopic 
galaxies from the SDSS; the second sample will be galaxies 
from the SDSS imaging catalogs. To the extent that grav- 
itational lensing can be neglected (justified in Appendix 
B), the only physical correlations arc at similar redshift, 
and so one can use the known redshift of one object of 
each pair to translate angles into transverse distances and 
fluxes into luminosities (Davis et al. 1978; Yee & Green 
1987; Phillipps & Shanks 1987; Lilje & Efstathiou 1988; 
Ferguson & Sandage 1991; Vader & Sandage 1991; Saun- 
ders et al. 1992; Lorrimer et al. 1994; Loveday 1997). For 
example, we will use the known redshift to select a fixed 
luminosity range for our imaging sample. 

It is well-known that the assumption of isotropic clus- 
tering permits one to deproject angular correlations into 
true spatial correlations (von Zeipel 1908; Fall & Tremaine 
1977; Davis et al. 1978; Phillipps 1985; Saunders et al. 
1992; Baugh & Efstathiou 1993; Loveday 1997; Dodelson 
& Gaztanaga 2000; Eisenstein & Zaldarriaga 2001). The 
method of Eisenstein (2003) performs this deprojection 
but notes that the noisy derivative that usually enters can 
be avoided by integrating the real-space spatial correlation 
function (denoted ^is) on spherical windows: 
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where W is our smoothing window and V = Airr^dr W{r). 
The resulting quantity A is physically useful: it is the av- 
erage overdensity of objects from the imaging catalog in 
the neighborhood (as defined by W) of a spectroscopic 
object. A of course depends on W{r), which in turn will 
have a characteristic scale a. Hence, A is a function of a, 
although wc suppress the label. 

Eisenstein (2003) demonstrates that the quantity A can 
be estimated simply by a pairwise summation over the two 
catalogs with a weighting function G{R), defined by 

GiR) = (2) 
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where the sums are over the spectroscopic and imaging 
catalogs, respectively. Nsp is the number of spectroscopic 
objects, Rjk is the transverse distance between the two 
objects (using the angle and the angular diameter distance 
to the given redshift), and (paizj) is the space density of 
the objects from the imaging catalog at the given redshift. 
This double summation can then be conveniently split, to 
define a value of A for each spectroscopic object: 

^j = lZ7^ E G{Rjk). (5) 



We can tabulate these A^ for each spectroscopic object 
and then take the average for any subset we wish. This 
allows us to find the run of overdensity versus redshift 
and luminosity and to perform jackknife resamplings with 
essentially no overhead. It also makes clear that there 
is no need to bin the spectroscopic galaxies in redshift 
before computing. Appendix A describes how we account 
for masks and efficiently handle the summation over large- 
separation pairs. 

Note that the space density of the imaging catalog as a 
function of redshift, (j)o{z), is not necessarily known. By 
choosing a sliding flux limit in the imaging catalog as a 
function of the spectroscopic object's redshift, we minimize 
the change in cf>Q (z) while also keeping the reference sample 
as constant as possible, so that clustering amplitudes can 
be most easily related to clustering bias. It is important to 
note that even if (poiz) is unknown, the relative values of 
A between different scales and different luminosities of the 
spectroscopic galaxies arc still accurately known. In other 
words, we can study differential effects at a single redshift 
and even the redshift evolution of those differential effects. 
Only the absolute clustering amplitudes as a function of 
redshift depend on the modeling of (j)Q{z). 

We will use windows of the form 



W{r) = exp 
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where a is a scale length. This choice allows us to focus 
on about one factor of two of scale in the cross-correlation 
^is(r). The apodization at small scale improves on a pure 
Gaussian window in that it downweights very small-angle 
pairs where object deblending might be a concern. 
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For this choice, the required F and G functions are 

F{R) = \c?se-^ [2s/o(s) - (1 + 2s)/i(s)] (7) 

G{R) = se-' [(6 - 8s)/o(s) + (-2 + 8s)7i(s)] (8) 

whore s = /Aa? and /q and Ii are modified Bessel func- 
tions of the first kind. For coding, it is useful to note that 

-1 
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for large s. We see that G « iR^ /2a^ at small R, so 
small separation pairs are indeed downweighted. For this 
window, F = 3a3(27r)3/2. 

For a power-law ^is(r) oc r~"^~", the resulting connection 
between A and is 
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For a = 0.8, this means that A = G.36^is(a) w 6s(1.76a). 
For a = 1.0, A = 6s(a)/3 « 6s(l-73a). In other words, 
we are roughly measuring the correlation function on the 
scale 1.75a. This shift is not surprising in that W{r) has 
most of its weight at r > o. 

We will compute the mean environment for many differ- 
ent luminosity bins, redshift bins, scale lengths, and two 
imaging galaxy samples. The errors on each data point 
arc estimated by jackknifo resampling. Wo use 50 spatially 
coherent subsamples of the data for the resamplings. We 
have investigated the covariances between different mea- 
surements. Of course, different scale lengths and imaging 
samples for the same primary LRG bins are highly covari- 
ant, typically 50-70% correlated between scales separated 
by a factor of two. However, we find that different lumi- 
nosity and redshift bins are generally close to independent. 
This is expected: there are not many LRGs in a given large 
non-linear structure. 

We will use proper (i.e. non-comoving) distances through- 
out this paper. The reason is that with this choice the 
computed quantity (^o^ is predicted to be redshift inde- 
pendent for the stable clustering ansatz (Peebles 1980). 
This is easy to see because (jioV/S. is the number of imaging 
galaxies surrounding the spectroscopic galaxy, and stable 
clustering predicts that these structures are time indepen- 
dent. On much larger scales, linear perturbation theory 
(e.g. Peebles 1980) predicts 00 A cx (1 + z)""!?^ for a fixed 
proper scale, where D{z) is the growth function (perhaps 
generalized to include the evolution of bias) and n is the 
power-spectrum spectral index P oc k""; for concordance 
cold dark matter cosmologies, this prediction is also close 
to constant if the bias is time- independent. Aside from the 
aesthetics of a constant baseline hypothesis, the choice of 
non-comoving distances also makes it simple to average 
over redshifts if the spectroscopic sample is imperfectly 
volume-limited. In practice, we average 4>[z)A rather than 
A, as the former is empirically close to constant in redshift. 

Appendix B argues that the effects of weak lensing are 
only of order 1% for the results in this paper. This is 
slightly below our best 1 a errors. We therefore neglect 
gravitational lensing. False correlations imprinted by pho- 
tometric calibration errors common to the two samples are 
shown Appendix A to be negligibly small for this applica- 
tion. 



3. SDSS SAMPLES 

3.1. Description of the SDSS 

The SDSS (York et al. 2000; Stoughton et al. 2002; 

Abazajian et al. 2003, 2004) is imaging 10"* square de- 
grees away from the Galactic Plane in 5 passbands, u, g, 
r, i, and z (Fukugita ot al. 1996; Gunn et al. 1998). Im- 
age processing (Lupton et al. 2001; Stoughton et al. 2002; 
Pier et al. 2003) and calibration (Hogg et al. 2001; Smith, 
Tucker et al. 2002) allow one to select galaxies, quasars, 
and stars for follow-up spectroscopy with twin fiber-fed 
double-spectrographs. The spectra cover 3800A to 9200A 
with a resolution of 1800. Targets are assigned to plug 
plates with a tiling algorithm that ensures nearly complete 
samples (Blanton et al. 2003a). 

Galaxy spectroscopic target selection proceeds by two 
algorithms. The primary sample (Strauss et al. 2002), 
referred to here as the MAIN sample, targets galaxies 
brighter than r < 17.77. The surface density of such galax- 
ies is about 90 per square degree. The LRG algorithm 
(Eisenstein et al. 2001) selects ~12 additional galaxies per 
square degree, using color-magnitude cuts in g, r, and i to 
select galaxies to r < 19.5 that are likely to be luminous 
early-types at redshifts up to ~ 0.5. The selection is ex- 
tremely efficient, and the redshift success rate is very high. 
A few galaxies (3 per square degree at z > 0.15) match- 
ing the rest-frame color and luminosity properties of the 
LRGs are extracted from the MAIN sample; we refer to 
this combined set as the LRG sample. In detail, there are 
two sections of the LRG algorithm, known as Cut I and 
Cut II and described in Eisenstein et al. (2001). 

We begin from a spectroscopic sample covering 3,836 
square degrees. The exact survey geometry is expressed in 
terms of spherical polygons and is known as lss_samplel4. 
This set contains 55,000 spectroscopic LRGs in the redshift 
range 0.15 < z < 0.55. 

3.2. The spectroscopic LRG sample 

The SDSS LRG sample is nearly volume- limited, but not 
perfectly so. In particular, the flux limit of the survey cre- 
ates a significant luminosity threshold at redshifts above 
0.37. We therefore focus on two volume-limited subsam- 
ples in this paper: 0.20 < z < 0.36 with -23.2 < Mg < 
-21.2, and 0.20 < z < 0.44 with -23.2 < Mg < -21.8. 
Here, the Mg is the rest-frame g-h&nd absolute magnitude 
at z = 0.3. This is computed from the observed r mag- 
nitude using the k and passive evolution corrections of 
the "non-star-forming" model presented in Appendix B of 
Eisenstein et al. (2001). We evolve the galaxies to z = 0.3 
rather than 2 = 0, in order to keep the results closer to 
the observations. See Zehavi et al. (2004c) for more details 
and plots of the number densities of these samples. The 
first sample totals 26,000 LRGs; the second totals 12,500 
LRGs, with about 6000 new ones. 

Note that because the angular cross-correlation method 
does not rely on comparisons of the positions of any two 
spectroscopic galaxies, the exact geometry or redshift se- 
lection function of the LRG sample is not required. How- 
ever, it is relevant to note that the LRG selection uses 
a luminosity cut that depends on rest-frame color. In- 
trinsically redder galaxies can enter the sample at lower 
luminosities (Eisenstein et al. 2001). This means that 
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our results at lower luminosities {Mg > —21.5) arc pref- 
erentially from the red edge of the red sequence of early 
types. However, essentially all LRGs are from the red se- 
quence, which is intrinsically quite narrow (Faber 1973; 
Visvanathan & Sandagc 1977; Bower ct al. 1992), indeed 
narrower than the photometric errors in a typical SDSS 
observation. Hogg et al. (2003a) computed the mean den- 
sity as a function of luminosity and color in the MAIN 
sample, where there is no color selection. They found no 
gradient of density with color at the high luminosity end. 
Hence, between the photometric errors across the red se- 
quence and the null result of Hogg et al. (2003a), we expect 
any gradient with respect to color to be small. 

In this paper, wc at times wish to express the Mg mag- 
nitudes in terms of luminosities relative to L* . We have 
adopted M* = —20.35 for this purpose. This number was 
derived from the value of M* = —20.44 in the °-^r band 
at = 0.1 (Blanton ct al. 2003b), applying an evolution of 
0.31 mag to z = 0.3, and converting from the "-^r band to 
the g band using a mid-type spectral energy distribution. 
We note, however, that this M* value is only approximate; 
the Mg magnitudes of the LRGs are the quantities closer 
to the observations. 



Table 1 
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NOTES. — "The r-band apparent magnitude m* we adopt as a func- 
tion of redshift to approximate a constant set of galaxies. The results 
roughly match the value of an L* galaxy at z = 0.3 in an Qm = 0.3 
flat universe. However, the model likely slides off of this reference, 
in the sense that at lower redshift we axe using more luminous (and 
hence less abundant) galaxies. 

^The tolerance in observed g — i to define a "red" galaxy. The color 
must be within 5(g — i) of the red sequence at the given redshift. 
''We use two samples, one from m* — 0.5 to ra* + 1.0, the other 
m* — 0.5 to m* + 0.4, denoted as the "1.0" and "0.4" samples, re- 
spectively. 

"^(p^z) is the comoving density for this magnitude range predicted by 
taking a representative sample from the SDSS MAIN sample (Blan- 
ton et al. 2002, 2003b) and moving the galaxies to the given redshift. 
Magnitude evolution of Q = 1.6 was assumed. The numbers in the 
table are in units of 10~^h^Mpc~^. Note, however, that elsewhere 
in this paper, we use proper densities and so (feo = {1 + z)^</i(z). 
''The fraction of those galaxies predicted to be classified as "red". 



3.3. Galaxies from SDSS Imaging 

Wc use 5305 square degrees of imaging from the SDSS 
to define our second sample of galaxies. The imaging com- 
pletely covers the spectroscopic sample, but the extended 
region is useful in that one can compute correlations to 
LRGs that arc near the spectroscopic boundaries (which 
are numerous because of the plate coverage of the SDSS). 
We extract all galaxies down to r = 21. 

At each redshift considered, we use only a fraction of 
galaxies. Ideally we would select a consistent set of galax- 
ies at all redshifts, but this is not possible in fine detail. 
Wc have attempted to select a fixed luminosity range, cor- 
rected for evolution. To do so, we define a reference galaxy 
and find its r band magnitude as a function of redshift in 
the usual cosmology, then select galaxies within a specified 
range relative to this magnitude. 

Our reference magnitude, denoted m* , is computed from 
an early-type galaxy spectral energy distribution, includ- 
ing k corrections and passive evolution. The resulting r- 
band reference magnitudes as a function of redshift are 
given in Table 1. We have scaled the model so that it 
matches our estimate of L* at z = 0.3. 

We consider two luminosity ranges, M* — 0.5 to M* -1-0.4 
and M* — 0.5 to M* + 1.0. The former allows us to use 
galaxies up to z = 0.44 at r < 21, while the latter reaches 
z = 0.36 (the approximate volume-limited redshift of the 
LRG sample) with more galaxies. We refer to these as the 
"0.4" and "1.0" samples, respectively. 

The m* values are indeed close to that of an L* galaxy, 
but not perfectly so. In particular, our early-type model 
does not evolve as much as recent luminosity functions 
(e.g., Q = 1.6, Blanton ct al. 2003b) and the red color 
makes m* fainter at higher redshift than a typical galaxy 
would. As a result, wc estimate that m* is actually di- 
verging from the actual track of an L* galaxy by about 0.2 
magnitudes per 0.2 in redshift, in the sense that at lower 
redshift we are using more luminous and less abundant 
galaxies. For typical luminosity functions in this luminos- 
ity range, every 0.1 mag of shift alters the densities by 
10%. Therefore, we expect that our samples may contain 
anomalous evolution at the level of (1 -|- zY'^. Given the 
mild trends of mean environment with luminosity around 
L* found in Hogg et al. (2003a), it seems unlikely that 
the modest luminosity shifts in our modeling would sig- 
nificantly alter the bias properties of the imaging sample 
across our redshift range. 

While we surely could construct a model that would 
match L* more carefully, it is not obvious what property of 
galaxies one should actually track. For example, since (j)oA 
is the number of selected galaxies near the LRG, one might 
be more interested in defining the sample so that cluster 
early-types were being consistently selected, which is likely 
to be closer to what we have done. Indeed, as galaxies 
near and far from LRGs may evolve differentially, there 
may be no truly comoving selection. Our view is instead 
empirical: Table 1 defines the apparent magnitude ranges 
that we have used, and redsliift-dependent interpretations 
should work to this definition. 

We have taken MAIN sample galaxies from lss_sainplel4 
to construct a luminosity-color function (Blanton et al. 
2002, 2003b) and used that sample to predict the lumi- 
nosity function of galaxies in the observed r band at each 
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redshift. The resulting comoving number densities are in 
Table 1. We assume luminosity evolution of Q = 1.6 (1.6 
magnitudes per unit redshift), independent of color, but 
we include no evolution of the spectral energy distributions 
themselves. As expected, the number density increases by 
about 20% per 0.2 in redshift. Most of this is due to the 
favorable k correction of blue galaxies boosting their num- 
ber in the r-band sample. 

We find that the SDSS luminosity fimction predicts a 
comoving number density of about (pcorn = O.Ol/i'' Mpc^'^ 
and (pcom = 0.045/i3 Mpc"^ for the "1.0" and "0.4" sam- 
ples, respectively. Note, however, we will be using physical 
distances for our a, so one should use (j>o = (/)com(l + z)^ 
when converting from (poA to A and hence to ^ig. 

Because LRGs arc often in clusters or rich groups, we 
will be interested in their correlations with red galaxies 
as well. We use the same r-band flirx cuts as above and 
then define a red galaxy as one that is within a particular 
tolerance in g — z color of the empirical color of the red 
sequence of early-type galaxies as a function of redshift 
in the SDSS. We set the tolerance in observed g — i to 
match an 0.08 mag blueward shift in rest-frame g — r color 
as a function of redshift; the resulting values are listed 
in Table 1. The considerable increase in this observed 
color range as a function of redshift is because the spectral 
energy distributions of red and not-so-red galaxies diverge 
as one moves blueward over the 4000A break. Note that 
our color cut is a mix of evolving and non-evolving models: 
the central color is the empirical one and hence is evolving 
with redshift, whereas the color range is derived from non- 
evolving templates. 

Table 1 lists the red fraction, i.e. the ratio of the num- 
ber density of the red galaxies to that of all galaxies in the 
r-band range, predicted by our modeling with the SDSS 
luminosity-color function. The red fraction is predicted 
to drop toward higher redshift despite the fact that there 
is no evolution of spectral energy distributions or differ- 
ential evolution of red and blue galaxies in our model- 
ing. The primary reason is that the fc-corrections of the 
blue galaxies are more favorable at higher redshift, such 
that for a given apparent r-band magnitude, one is looking 
further down the the luminosity function of blue galaxies 
than that of red galaxies. A secondary reason is that in 
this model galaxies are evolving in luminosity somewhat 
faster than our m* reference value, so that one is selecting 
slightly lower Imninosity and hence slightly bluer galaxies 
at higher redshift. 

Despite our effort to match rest-frame colors, it is possi- 
ble that our selection of red galaxies is redshift dependent. 
Modeling errors in some aspect of the color selection or in 
the luminosity evolution of the galaxy population would 
be one cause; another would be the neglect of photometric 
scatter across the color boundaries. As before, our sample 
definition is empirically precise, but modelers will need to 
account for our selection. 

To conclude, we stress that while uncertainties in the 
evolution of our samples must be addressed to interpret 
bulk redshift trends, the uncertain 4>o values cancel out 
and are no concern when studying trends with scale and 
LRG luminosity within volume- limited samples such as are 
employed here. In addition, one can study the redshift 
evolution of trends with scale and luminosity. 




Fig. 1. — Mean overdensity {(f>oA) around red galaxies as a func- 
tion of luminosity, expressed as the rest-frame g magnitude, k+e cor- 
rected to z = 0.3. The window is o = 0.125ft~^ Mpc, proper. Only 
galaxies in the range 0.20 < z < 0.36 have been used. The "1.0" 
imaging sample is used, so the galaxies are in the range M* — 0.5 
to M* + 1.0. 



4. RESULTS 
4.1. Scale and Luminosity Dependences 

In Figure 1, we present (fioA with a proper scale length of 
a = 0.125ft,-iMpc for LRGs in the redshift range 0.20 < 
z < 0.36 with regard to imaging galaxies in the range 
M* - 0.5 to M* + 1.0 (the "1.0" sample). The LRGs have 
been binned into 0.1 magnitude bins from Mg = —21.2 to 
-22.7 (i.e. 2.2L* to 8.3L*). It is immediately clear that the 
density of galaxies around the LRGs is a strong function 
of LRG luminosity, ranging from ^qA = 8 to ^30. 

What do these numbers mean? (pQ is the proper number 
density of the imaging sample galaxies, ~ 0.022/i^ Mpc~^ 
at z = 0.3. Hence, we have A running from ^ 400 to 
~ 1500. This is the value of the cross-correlation function 
between LRGs and ~ L* galaxies, averaged over the win- 
dow Wir) and centered at ~0.22ft,-i Mpc. Note that the 
classical £,{r) = (r/ro)^^'^ with rg = 5/i^^Mpc comoving 
would predict A = 170, but LRGs are known to be highly 
biased (White, TuUy, & Davis 1988; Norberg et al. 2001; 
Zehavi et al. 2004b,c). 

Alternatively, if one multiplies 4>oA by the volume of 
W{r) (here, V = 3(27ra2)i-5 = 0.0922/1^3 Mpc^), ^-^^^ 
has the average number of M* — 0.5 to M* + 1.0 galaxies 
(above a random unclustered floor of (poV, insignificant for 
a < 0.5/i~^Mpc but not for a > l/i~^Mpc) surrounding 
the LRG, with the counting weighted by W{r). This latter 
interpretation is independent of (po and so our numbers 
are quite precise, e.g. for = 10.8, there is 1.0 galaxy, 
weighting the count by W{r), in the range M* — 0.5 to 
M* + 1.0 around that class of LRGs. Note that W{r) < 1, 
so the weighted count is less than the actual number of 
galaxies. We will show 4>oA when we want to appeal to 
the correlation function and (poVA when we want to stress 
the empirical accuracy. 

More approximately, if one considered the density of 
galaxies to be constant in the window region, then cp^A 
would be that density (proper, in this case). In other 
words, the proper density of M* — 0.5 to M* + 1.0 galaxies 
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0.2/1-1 Mpc from an LRG is about lO-SO/i^ Mpc'^ Of 
course, since the density is in fact steeply declining with 
scale, this interpretation is not precise. 

One might worry that because (f)o is a reasonably strong 
function of z (sec Table 1 and recall the extra factor of 
(1 + z)^), then deviations in the mean redshift as a func- 
tion of luminosity would compromise the comparison of A 
from Figure 1. In fact, the samples are volume- limited to 
sufficient accuracy: the mean redshift of the higher lumi- 
nosity bins in Figure 1 is only 0.004 higher than that of the 
lower luminosity bins (0.296 vs. 0.292), which is only 1.2% 
more in (1 -|- z)^. This demonstrates that cpo is essentially 
constant across the LRG luminosity range being used. 

With the result that (po is the same for all of the lu- 
minosity bins, Figure 1 shows a four-fold variation in A 
as a function of LRG luminosity, which means a four-fold 
variation in the small-scale cross-correlation of LRGs with 
respect to L* galaxies! 

It is difficult to quote an exact value of the bias of the 
LRGs with respect to the mass because of the uncertainty 
in ^0 and the fact that we are using a particular set of 
galaxies (namely, L* galaxies, as defined by Table 1) to 
trace the density field. Formally, we arc measuring the 
cross-correlation between LRGs and L* galaxies. There is 
the temptation to interpret this cross-correlation as vary- 
ing as the product of the bias of LRGs with respect to mass 
and the bias of L* galaxies with respect to mass. Since the 
latter is constant, one would interpret the change in A as 
only the change in the bias of LRGs. However, galaxy bias 
need not be separable in this simple fashion. Although L* 
galaxies seem to have bias close to unity whole (e.g. 
Verde et al. 2002), the relation of their numbers to mass 
might be non- linear, yielding a different slope in the more 
extreme environments traced by the LRGs than the slope 
found when averaging over all environments. 

Figure 1 is not quite a linear relation. If we instead 
plot the density versus the luminosity raised to the 1.5 
power, then we find a tight linear relation with a non- 
zero intercept. This is shown in the first panel of Figure 
2. The horizontal axis reads as luminosity relative to L* 
(M* = —20.35), but we have stretched the axis into a L^'^ 
dependence. We find that the fit AL^-^ -|- i? is a good fit in 
all cases, and it is slightly better (2-a) in most cases than 
AL + B, particularly at larger scales. The per degree 
of freedom is typically about 1, taking the jackknife errors 
in different luminosity bins to be independent. We will 
therefore quote our quantitative results as the parameters 
of this fit. 

Figure 2 shows (j)oVA as a function of LRG luminos- 
ity for 6 different scales, increasing by factors of two from 
a = 0.125ft-iMpc to a = 4:h~^ Mpc, proper. We have 
overplotted the best AL^-^+B fit in each case. The param- 
eters of the best-fit line for these and other subsamples are 
given in Table 2 and 3. We quote the intercept of the lin- 
ear fit at an intermediate luminosity value chosen so that 
the error in the intercept is minimized and is not covariant 
with the estimate of the slope. The errors on the slopes 
indicate that the luminosity dependence of clustering is de- 
tected at about 20-tT over the range —22.7 < Mg < —21.2. 

Figure 2 and Table 2 show that there is a significant 
change in slope as a function of scale. This is shown 
directly in Figure 2d, where the best-fit lines from a = 



0.125/1-1 Mpc and a = Ah^^Mpc have been rescaled and 
overplotted. Larger scales have a softer run of A against 
luminosity; for example, a = Ah~^ Mpc, the overdensity 
varies only by a factor of 2 across the same range of lumi- 
nosity. Turning this around, this means that high luminos- 
ity LRGs have a steeper scaling of A with scale. In Table 
2, we report a slope of 0.110 ± 0.005 for a = 0.125/i-i Mpc 
and 0.058 ± 0.010 for a = 2h-^ Mpc; treating these two 
widely separated scales as independent, this means that 
the scale dependence of the luminosity dependence is de- 
tected at 4.6-(T. 

We next consider the detailed shape as a function of 
scale a. Figure 3 shows the value of cf){)lS. as a function 
of a for 3 coarser bins in luminosity. The bottom panel 
shows the residuals relative to an power law. Devi- 
ations from the power-law model are obvious. Using the 
full covariance matrices, the best-fit power laws have, from 
high to low luminosity, = 18, 56, and 88 for 4 degrees 
of freedom, and hence power laws are strongly rejected 
with goodness-of-fit probabilities of 10"'^, 2 x IQ-^i, and 
3 X IQ-i®, respectively. The deficit at ^ lh~^ Mpc and ex- 
cesses at smaller and larger scales match the behavior seen 
in the SDSS MAIN sample in Zehavi et al. (2004a) and in 
the LRG sample in Zehavi et al. (2004c). It is interest- 
ing to interpret these correlations in the halo occupation 
model as a sum of two terms, a small-scale component in 
which both galaxies are in the same halo and a large-scale 
component in which both galaxies are in different halos. 
Clearly, both terms would be required to fit Figure 3. 

The fact that the high and low luminosity curves in 
Figure 3 are closer together at larger scale than at lower 
scale is another manifestation of the mild scale dependence 
of the Imninosity dependence. 

As stated in §2, measurements at different scales are 
covariant. We have used our 50 jackknife samples to con- 
struct the 6-dimensional covariance matrix for Figure 3 
and its "0.4" imaging sample equivalent. With only 50 
jackknife samples, the results are noisy, but the covari- 
ance between scales appears to reduce from unity by about 
50-70% per factor of two in scale (i.e., nearest scales are 
50-70% covariant, next nearest are 25-49% covariant, etc.). 
The lower-luminosity samples, i.e. — 21.7 < Mg < —21.2, 
are more covariant (70%), whereas the more luminous 
LRGs are less covariant (50%). Widely separated scales, 
such as in the 0.125/1-^ Mpc to 2h~^ Mpc comparison above, 
are at most 20% covariant. 

Tables 2 and 3 show mild evidence for a 10% decrease 
of ^oA with increasing redshift (in the non-color-selected 
samples). We see no significant evidence for a change in 
slope of A versus L^^^ as a function of redshift. As stated 
in §2, both stable clustering and linear theory predict a 
nearly redshift independent result. Unfortunately, uncer- 
tainties in the redshift evolution of the imaging sample 
currently prevent us from describing the exact evolution 
of (/'qA and confirming these hypotheses. If the results 
from the extrapolation of the 2; = 0.1 SDSS luminosity- 
color function (Table 1) are correct, then the comoving 
density of the imaging sample is increasing slightly with 
redshift, such that the 4>oA for a sample of fixed comov- 
ing density would be scaling as {1 + z)~^. This would 
argue for some anomaly beyond stable clustering, perhaps 
caused by some kind of unmodeled evolution. However, it 
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Fig. 2. — Mean ovordcnsity ((j>oVA) around red galaxies as a function of luminosity, expressed relative to M* = —21.35. The horizontal 
axis has been stretched to scale as luminosity to the 1.5 power, as this makes the linear fit better in most cases. The solid lines are the 
best AL^'^ + B fit, reported for these and other samples in Tables 2 and 3. (a-f) Six panels show different scales; a = 0.125/i~^ Mpc, 
a = 0.25h~^ Mpc, a = 0.5h~^ Mpc, a = lh~^ Mpc, a = 2h~^ Mpc, and a = 4h~^ Mpc, proper, respectively. The dashed lines in panel (d) are 
the best-fit lines from a = 0.125fe~^ Mpc (steeper) and a = 4h~^ Mpc (shallower) rescaled and overplotted to show how the slope is changing 
with scale. Recall that the A statistic is in fact probing the correlation function on scales of about 1.75a, so the range here is roughly 0.2h~ 
to 7/i-i Mpc. Only galaxies in the range 0.20 <z< 0.36 have been used. The "1.0" imaging sample is used (M* - 0.5 to M* + 1.0). 
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Fig. 3. — Mean overdensity around red galaxies as a function 
of proper scale in three different luminosity bins. The top panel 
shows the quantity cpoA; the bottom panel shows a'^(f>oA, which is 
chosen to flatten out the curves so that one can see the fine detail. 
The deviations from a power-law fis are clearly visible, as is the 
luminosity dependence to the clustering. More subtle is the fact 
that the luminosity dependence is slightly stronger at the smaller 
scales. The reader is reminded that the A statistic measures the 
correlation function on scales of ~ 1.75a. Only galaxies in the range 
0.20 < z < 0.36 have been used. The "1.0" imaging sample is used 
(M* - 0.5 to M* + 1.0). 

may be that cluster galaxies evolve similarly to our model 
and that it is the field galaxies that are driving the evo- 
lution of the luminosity function (e.g. Lin et al. 1999), in 
which case stable clustering would predict constant (j)oA. 
regardless of the bulk evolution in (poi^)- Calibrating our 
empirical measurement will require more precise model- 
ing of the redshift distribution of SDSS galaxies at these 
flux levels. In principle, we could measure the differential 
evolution between different scales, but the predicted varia- 
tions are smaller than our measurement uncertainties. On 
the positive side, the insensitivity to redshift validates our 
averaging (j)oA over galaxies at many redshifts; any devia- 
tions from a volume-limited LRG sample (which are minor 
in any case) would produce negligible biases in A. 

4.2. Fraction of red galaxies 

We define the red fraction as the ratio of (poA for our red 
galaxy imaging sample to that for the full imaging sample. 
Because (poVA+ipoV really is the average number of galax- 
ies, weighted by W{r), surrounding our LRGs, this ratio 
of these quantities is the fraction of those galaxies that are 
red, as defined in §3.3. We neglect the homogeneous term 
(poV because it is negligible at small scales (only 10% of 
the clustered term at a = Ih^^ Mpc and dropping as 
below that), whereas the red fractions at larger scales are 
close enough to the field values that including the homo- 
geneous term wouldn't change the ratio appreciably (e.g., 
for a = 4:h~^ Mpc, the results might drop by 3%). 

Figure 4 shows the red galaxy fraction as a function 



of scale and LRG luminosity. The red galaxy fraction is 
a strong function of scale a, as one would expect from 
the excess of early-type galaxies in dense regions (Abell 
1965; Oemler 1974; Melnick & Sargent 1977; Dressier 1980; 
Postman & Geller 1984). We see a slight increase in the 
red fraction around higher luminosity LRGs, but only on 
small scales. On scales above 1 Mpc, we do not detect a 
trend with luminosity. The higher luminosity "0.4" imag- 
ing sample has a slightly higher red fraction, as one would 
expect from the color-magnitude distribution of galaxies. 

We display two different redshift bins in the bottom pan- 
els of Figure 4. The higher redshift bin has a smaller red 
fraction. While this trend is in the same direction as the 
conventional wisdom that galaxies are bluer at higher red- 
shift, it is also possible that our selection definition has 
created a moving target, as discussed in §3.3 One effect 
that is certainly present is that our r-band fiux limits 
select galaxies at bluer rest-frame wavelengths at higher 
redshift, causing the samples to tilt toward a bluer frac- 
tion. In addition, our definition of a "red" galaxy could be 
redshift dependent due to imperfect modeling of the rest- 
frame colors. It is interesting to note that the red frac- 
tions at the a = 4/i~^ Mpc scale are only slightly higher 
than the fractions predicted from the low-redshift SDSS 
luminosity-color function in Table 1. This suggests that at 

Ih"^ Mpc (10/i^^Mpc comoving), one has nearly con- 
verged back to the field value, despite having chosen the 
galaxies by their proximity to a LRG. 

The top panel of Figure 5 shows the red fraction as a 
function of scale for 3 redshift bins and the "0.4" imaging 
sample. Again, the red fraction is redshift dependent. In 
the bottom panel, we attempt to correct for the simplest 
effect, namely that lower luminosity blue galaxies can en- 
ter the sample at higher redshift simply because of more 
favorable fc-corrections. We do this by simply reducing 
the number of blue galaxies at higher redshift by some 
fraction x {x = 1 being no change). If the red fraction 
at the higher redshift is /2, then the fraction /i at the 
lower redshift, having diluted the blue galaxy density by 
X, will be /i = {1 — X + x/ f2)^^- Note that this is not 
simply a multiplicative scaling in the red fraction. From 
the model in §3.3 and Table 1 in which we extrapolate the 
low-redshift luminosity-color function without any color 
evolution, we derive x = 0.86 (0.75) for 0.28 < z < 0.36 
(0.36 < z < 0.44) relative to the 0.20 < z < 0.28 sample. 
Applying these corrections, we see that the red fractions 
at large scales (i.e., the field) are consistent with being in- 
dependent of redshift. On smaller scales, there is a trend 
with redshift, suggesting that in regions near LRGs (i.e., 
high-density regions), galaxies have reddened with time. 
Strictly speaking, this is the evolution of a population 
defined by location rather than one that is consistently 
tagged across time. 

We stress again to the reader that variations in the def- 
inition of "red" as a function of redshift could still move 
the curves up or down, so one should not conclude that 
there is no evolution of the red fraction in the field. These 
uncertainties should be amenable to better modeling of the 
galaxy fc-corrections and luminosity functions. For exam- 
ple, differential luminosity evolution between red and blue 
galaxies (Lin et al. 1999) could alter the red fraction pre- 
dictions from Table 1. Nevertheless, the conclusion that 
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Fig. 4. — The fraction of red galaxies in the imaging sample as a function of LRG luminosity for different scales. Each line is a different 
scale, from a = 0.125fe~^ Mpc to 4fe~^ Mpc by factor of two, from top to bottom. A red galaxy is defined as being within 5{g — i) mag in 
observed g — i color from the red sequence color at the LRG redshift (see Table 1); this classification is not redshift independent and so the 
fraction of red galaxies is hard to quantify versus redshift. However, because the samples are volume-limited, comparisons between scales 
and LRG luminosities are well posed. The errors arc computed from the variations in the red fraction among the jackknifc samples, not from 
propagation of errors in (/iqA. (Upper left) M* -|- 1.0 sample for 0.20 < z < 0.36. (Upper right) M* + 0.4 sample for 0.20 < z < 0.36. (Lower 
left) M* + 1.0 sample for 0.20 <z < 0.28. (Lower right) M* + 0.4 sample for 0.28 <z< 0.36. 
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Fig. 5. — {Top panel) The fraction of red galaxies in the imaging 
sample as a function of scale for different redshift bins. In all cases, 
only LRGs in the range —22.2 < Mg < —21.8 are used, and the 
imaging galaxies are from the "0.4" sample. The points with error 
bars are the measurements with 0.20 < z < 0.28, 0.28 < z < 0.36, 
and 0.36 < z < 0.44, from top to bottom. {Bottom panel) The same, 
but the two higher redshift bins have been rescaled according to the 
model that the blue galaxies are to be diluted by 86% and 75%, 
respectively. These numbers were determined from extrapolation of 
the MAIN sample in the absence of any evolution of intrinsic color. 
The red fractions on larger scales are consistent after this rescaling; 
the smaller scales show a significant difference. 



the evolution of the red fraction is different between re- 
gions near LRGs and far from LRGs is robust, as we have 
used the same definition of "red" at aU scales. 

We also caution against interpreting galaxies near LRGs 
as necessarily representing a cluster (e.g., mass above lO^^M©) 
population, as not all LRGs live in clusters (e.g., Loh et 
al. 2003). However, the densities in these regions are high 
on average, about le/i^Mpc"^ for the a = 0.125/1"^ Mpc 
case, and so our results will surely bear on the question of 
differential evolution between high and low density envi- 
ronments (e.g. Hashimoto et al. 1998; Balogh et al. 1999; 
Margoniner et al. 2001; Lewis et al. 2002; Gomez et al. 
2003; de Propris et al. 2004). 

5. DISCUSSION 

We have demonstrated that the environments of LRGs, 
as measured by the surrounding overdensity of L* galax- 
ies, varies strongly with luminosity. Across the range from 
2.2L* to 8.3L*, the clustering amplitude changes by a fac- 
tor of 4 on 0.2h^^ Mpc scales and a factor of 2 on 7h~^ Mpc 
scales. This trend was clearly seen before (Norberg et al. 
2001; Zehavi et al. 2002; Hogg et al. 2003a), but here we 
use the large volume of the LRG sample to bring this result 
to very high signal-to-noisc ratio. Moreover, the variation 
of the slope versus luminosity as a function of scale im- 
plies that LRGs have a scale-dependent bias that varies 
with luminosity. Higher luminosity LRGs have even more 
clustering on sub-Mpc scales than one would project from 



their clustering on ^ 5h~^ Mpc scales. 

The variation of clustering amplitude as a function of 
LRG luminosity is striking in its strength, and of course 
these galaxies are already more clustered than less lumi- 
nous galaxies. This clearly points to a significant change 
in the masses of the host halos of LRGs as a function of 
luminosity. 

The clustering as a function of scale shows clear devia- 
tions from a power law, in a manner quite consistent with 
the LRG autocorrelation function (Zehavi et al. 2004c) 
and with less luminous samples (Zehavi et al. 2004a, b). 
The interpretation in terms of the one-halo and two-halo 
terms of the halo occupation model (Ma & Fry 2000; Pea- 
cock & Smith 2000; Seljak 2000; Scoccimarro et al. 2001; 
Berlind & Weinberg 2002) would naturally explain the 
qualitative structure (Berlind et al. 2003; Magliocchetti 
& Porciani 2003; Scranton 2003; Zehavi et al. 2004a); we 
will pursue quantitative analyses in future papers. Com- 
bining the cross-correlations methods employed here with 
the LRG and L* galaxy auto-correlations should yield new 
constraints on the LRG and L* galaxy populations in mas- 
sive halos. 

The fraction of red galaxies as a function of scale and 
luminosity are qualitatively consistent with the familiar 
density-morphology relation (Dressier 1980; Postman & 
Geller 1984). Higher luminosity LRGs are surrounded by 
a slightly larger fraction of red galaxies at sub-Mpc scales. 

Tracking redshift evolution is a challenge to any cluster- 
ing method, as it requires detailed understanding of the 
evolution of the selection. In our case, this is phrased as 
the need for a measurement of the number density 0o of the 
imaging sample. Redshift evolution of differential effects 
can still be measured robustly. We find, for example, that 
the relation of red fraction versus scale changes with red- 
shift, such that the high-density regions near LRGs have 
reddened more than low-density regions far from LRGs. 
This result cannot be avoided simply by rescaling the num- 
ber of blue galaxies, thereby addressing the most obvious 
gap in our modeling. We intend to pursue this issue fur- 
ther in future work. More generally, with precise modeling 
of the redshift distribution of the galaxies in the SDSS im- 
ages, we could recover 0o and measure the evolution of 
A. 

We propose that the scale-dependent clustering of LRGs 
could serve as a test of the formation theories for these 
galaxies. In particular, the strong clustering on small 
scales is a challenge for theories that rely on passive evo- 
lution from high redshift, with no environment-dependent 
evolution of the galaxy. If LRGs sit in only the most mas- 
sive halos, then they will be highly clustered (Kaiser 1984; 
Bardeen et al. 1986; Mo & White 1996; Sheth, Mo, & Tor- 
men 2001), but the correlation between their halos' mass 
today and their environment at high redshift may not be 
sufficiently tight to allow the high-redshift environment 
to dictate today's luminosity. Processes that enhance the 
LRG luminosity in situ would loosen these constraints. 

More generally, this work highlights the considerable ad- 
vantage of using cross-correlations between imaging and 
spectroscopic data sets for the study of the small-scale 
clustering of rare types of galaxies. On small scales, one 
is nearly always shot-noise limited, and auto-correlation 
analyses will suffer two portions of shot noise for rare ob- 
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jects. Cross-correlations with more populous sets of galax- 
ies only incur the shot noise once. Of course, one would 
expect that the clustering trends arc reduced in the cross- 
correlation, but this is a mild loss of signal compared to the 
gains in the signal-to-noisc ratio. Moreover, for luminous 
galaxies, acquiring spectroscopic redshifts for the fainter 
L* companions is doubly challenging, as they are both 
faint and numerous. Angular cross-correlations therefore 
allow one to leverage a considerable amount of imaging 
data. Bcc;ause the correlations cfFectivcly impose a redshift 
on the imaging data, one can extract rest-frame properties 
from the imaging set (such as the red fractions presented 
here) despite not having true redshifts. 

We expect that this method could have considerable ap- 
plication to deep wide-field data from the next generation 
of ground-based surveys or from the Spitzer satellite, as it 
offers a precise and quantitative measurement of clustering 
with a minimum of spectroscopic information. 
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APPENDIX A 

IMPLEMENTATION DETAILS 

When implementing the summation for the Aj, it is 
easy to take account of masks and boundaries as well as 
to avoid summing over very widely separated pairs. Fol- 
lowing Eisenstein (2003), we limit the explicit summation 
at some Rmax (we use 10a) and add the remaining piece 

— /"OO 

Pouter = / ^-KRdR [1 + wM\G{R) (Al) 

Here, n is the number density of objects from the imaging 
catalog (selected at the chosen redshift) per unit area on 
the sky (measured in transverse distance units). The ho- 
mogeneous term (the 1 in the square brackets) has an easy 
integral involving F{R), whereas the correlated term can 
be easily done using a power-law approximation for 
and the large-radius expansion for G{R). If oc r~^~" 
so that Wis oc r~", then we have 



^Wis{R) = R^UR)^^^^^ 



(A2) 



and so 

^outer 



n 2ttF{R^^^) 
(^n V 



■ ^is(-^max) 



r(f)v^ 



r(«±l)^(l + a)- 



(A3) 

The second term depends on ^is, which is what we are 
trying to find via A. However, if the correction is small, 
then it is a good approximation to keep the shape of 
fixed while trying to find its amplitude. In this case, the 
second term becomes 

A _ ^ 6s(-Rmax) P I \A\ 

^outer,clust — — ^ _^-| , (A4j 



where 



0- 



r(f)V^3^/7^. 



^is(a) 7r(l-|-a) 
1 



(A5) 



T{^) 2 V 2 ' " r(2 - f) 

f3 = 10.49 for a = 0.75 and /3 = 37r for a = 1. 

Hence, to include the effects of the upper limit to the 
explicit summation, we now have 



A,- 



1 



iV 



-2nnF{R 



Cis (^n 



max; ~t~ ^ ] 

ke{im},R.ik<R 
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G{Rjk) 



~ Ua) Hl + a) ^^^^ 

We now would average the Aj over many spectroscopic 

objects to estimate A. However, note that the last term 
is just linear in A, so we can instead compute 



A', 



1 



(poV 



max 

ke{im},Rjk<IiaiBx 



for each object, and then estimate A for any subset as 



(A7) 



A 



1 



(A8) 



{i=(a) 7r(l-|-a) 

where the angle brackets indicate averaging over a subset 
of spectroscopic objects. 

An inner integration limit can be treated the same way, 
although the integral for the correlated term is different. 

dinner = T"?? / 2wRdR [1 + Wis{R)]G{R) (A9) 



n 2nF{R„ 



r 

Jo 



V 



RdRG{R) 



^ 27ri?^i„Cis(i?min) r(f )V^ 



V 



r(^) 



/O V 

For the W{r) in equation 6, the second term becomes 

5 



(AlO) 



Hnner,clust 



= A 



R„ 



a 



Cis (-^min) 



/3 



6s (a) 



(All) 



/87r(4 - a) 

Because of the apodization of the window at small radii, 
this term is quite small, typically of order {Rmin/a)^- How- 
ever, other choices of W{r) that have W{0) ^ will have 
larger contributions. 

Finally, the same principles apply to masked regions. 
As described in Eisenstein (2003), it is easy to solve the 
problem by Monte Carlo by creating a dense set of ran- 
dom points on the sky outside of the survey (i.e., filling 
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the masks and any border regions). We then add to each 
Aj a quantity computed by summing over the random 
points found within the Rmm to i?max radial interval, sub- 
ject to the weighting G{R)[1 +Wis{R)]{n/nrandom), where 
firandom IS the arcal density of the random catalog in the 
same units of n. Again, the correlated term can be fi- 
nessed by approximating Wis{R) = Wis{a){R/a)~°' , as it 
then becomes 



^mask.clust 



Vn. 



random 



-A/3 



E 



R 



Hence, in practice, we compute 



(A12) 



A' 



+2TTn{F{R^in) - -F(i?max)) 



+ 



n 



"^random 



E 



kZ{ran},Rjk <Rt, 



and 



a(3 



G{Rjk) 



E GiR,,) ( ^ 



(A13) 



^l^random , ^ , . „ \ CI 



(A14) 

For a particular subset of spectroscopic objects, we esti- 
mate A as 



Ci,(a) 



JL 



/87r(4-a) 



(^) 



5 ^^(flm.n) 



(A15) 

Again, /3 « 10 and a is typically 0.7 to 1. Obviously, one 
would like to pick i?max and Rmin so that the corrections in 
the denominator are minimal; we use i?max = 10a, which 
should be less than 2% correction, and i?min = 30/i~^ kpc, 
which is a tiny correction. We exclude any objects for 
which the mask covered more than 40% of the iimax cir- 
cular region, but in practice most objects have very lit- 
tle masking. We find (Mj) to be 1% or less, even for 
a = 4:h~^ Mpc with i^^ax = 40/i~^Mpc. For i?max = 
40h~^ Mpc, about 1/3 of the objects have some fraction of 
their region masked; this number drops rapidly for smaller 

-^max* 

Our formulae use the flat-sky approximation, which is 
well-satisfied since the angle on the sky even for the largest 
-Rmax (40/i~^Mpc) and closest LRGs is only 5° and the 
angles contributing significantly to the integrals are yet 
smaller. The residuals would be of order 6"^, which is 
<C 1%. We note that it is important to define the fiat- 
sky radius R as 2D sm{9/2) for distance D to the LRG 
and separation angle 9 between the two objects, as this 
preserves a homogeneous distribution from the sphere to 
the flat sky and permits the unclustered background to 
cancel exactly from A. 

At present, our mask only includes the major boundaries 
of the survey, not small-scale effects such as bright stars or 
bad columns. If these masks were uncorrelated with the 
catalog of spectroscopic objects, then the neglect of these 
regions would not bias the values of A. Of course, the 
masks and catalogs are correlated. We assess the amount 
of spurious correlation by using an catalog of M stars as 



our "spectroscopic" set. We choose reasonably bright, ex- 
tremely red M stars, as this region of color-magnitude 
space has essentially no galaxy contaminants that could 
correlate with the imaging catalog. We place the stars fic- 
titiously at z = 0.3 and compute the correlations with the 
galaxies. The results are consistent with zero, with errors 
that are roughly the same as the errors in the LRG clus- 
tering signal. In other words, the neglect of the small-scale 
mask is less than a 1-a effect for our measurements. 

Photometric calibration errors would correlate the spec- 
troscopic LRGs with the imaging catalog, since both sam- 
ples impose a flux limit, and create a false signal. However, 
the effect is small. Following equation (Al), we have 



yAcai = n 



2t^R dR Wi^G{R), 



(A16) 



where Wis is now the cross-correlation between the two 
samples due to the calibration error. If Wis were scale- 
independent for scales near a, as would be the case if the 
calibration were wrong in patches whose size was much 
larger than o, then this integral would cancel to zero. To 
avoid cancellation, one must have structure in Wis near 
the scale a. The hmction G{R) peaks at 0.8 &i R k, a. 
Conservatively, one would have 0o^Acai = n'Ka?Wis{a). 
n « 8/i^Mpc"^ for the "1.0" sample aX z = 0.3, so for 
a = l/i"-'^ Mpc, we have (/)oV^Acai ~ 25i«is. The observed 
^o^A exceeds 5 with 3% errors. Hence, we need the 
calibration-induced Wis to be less than 0.006. The SDSS 
rms error in the r band are about 2%, which produces a 
6% response in the LRGs (Eisenstein et al. 2001) and a 2% 
response in the imaging sample. This implies = 0.0012, 
a factor of 5 smaller than required. This is conservative 
because the actual correlation function from the SDSS cal- 
ibration errors is a smooth function of scale, i.e., the 2% er- 
rors accrue from errors on many patch sizes, which tends to 
cause the integral in equation (A16) to be smaller. Hence, 
angular correlations between the samples due to calibra- 
tion errors are negligibly small for our purpose. 

It is worth noting that the deprojection formalism can 
be applied trivially to spectroscopic auto-correlation ap- 
plications with the usual projected correlation function 



w, 



/oo 
dZ ^rediVZ^ + R^), (A17) 
-oo 

where ^redi^) is the redshift space correlation function. 
Defining A as in Equation (1), we have 

A=^ f 2nRdRwp{R)G{R) (A18) 

y Jo 

Hence, to compute A from an auto-correlation study, there 
is no need to deproject Wp to ^reai and then integrate to 
A. 



APPENDIX B 

WEAK LENSING EFFECTS 

We argue here that weak lensing magnification effects 
are small for our analysis, only a ^ 1% effect, compa- 
rable to our quoted errors in the best cases. The root 
reasons are simple: our flux limits select L* galaxies at 
the redshift of the LRG, so the luminosity implied for a 
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target at a redshift high enough to be lensed effectively 
is much higher, resulting in a low density of potential 
sources. In more observational terms, at ^; = 0.3 we are 
using r = 20 galaxies and there are rather few r = 20 
galaxies at z > 0.6 on the sky. Furthermore, the scales we 
are studying (> 0.2ft,~^ Mpc) are large compared to the 
Einstein radius of an LRG or even a reasonable cluster. 

Mathematically, we consider magnification patterns of 
the form {R/Re)~^, where Re is a constant. This would 
generate an angular correlation proportional to i?^^, with 
the amplitude depending on the cosmological distances 
and the slope of the luminosity function. 

For example, summing over source redshifts, one could 
write the contribution to A as 

^ ^ 2TTRdRG{R) I dz(j){z) 

J ZLRG 

2 



WL 



V 



Dm{z) 



{a -I) 



Re{z) 
R 



(Bl) 



_Dm{zlrg) . 

where z is the distance along the line of sight (not redshift), 
R is the transverse distance across the line of sight, Dm{z) 
is the cosmological proper motion distance, (f){z) is the 
comoving density of galaxies passing the selection at the 
LRG redshift, a is the logarithmic slope of the cumulative 
luminosity function, and Re{z) is the Einstein radius (in 
distance, not angle) for that source redshift. The ratio of 
Dm enters because of the volume per solid angle, and the 
a — 1 is the familiar competition between magnification of 
sources and dilution of source density (Turner ct al. 1984). 

The equation can be simplified by noting that equations 
(10), (A2), and 



A = 



2'KRdRwis{R)G{R) 



[see eq. (Al)] imply that 
1 
V 



2TrRdR G{R) 



Ro 



Ro 



(B2) 



(B3) 



for the W{r) in equation (6). This means that the weak 
lensing contribution to A is 



1 



37ra^ 



f 

J ZT 



dz 



Dm{z) 
Dm{zlrg) 



{a-l)RE{z). 

(B4) 

One could do a careful integration using a luminos- 
ity function, but here we will merely estimate. Re is 
small unless the sources are well behind the lens. Char- 
acteristically, we would have z ^ 2zlrg ^md Dm{z) ~ 
2Dm{zlrg)- The luminosity threshold will rise as the 
square of the luminosity distance, which is (1 + z)Dm, 
in addition to any k corrections, so we would typically 
have the luminosity threshold increased by about a fac- 
tor of 5 10. This suggests (p{z) w 0.01(f)o- We assume 
Q — 1 ~ 1. The line of sight interval would characteris- 
tically be 1000/i~^Mpc. We conservatively take Re = 
20 kpc, which is already cluster sized (6" radius) and 
would correspond to 10% magnification ratios at 1' (e.g., 
Bartelmarm & Schneider 2001; Benitez et al. 2001; Jain 
et al. 2003). Putting these together suggests Awl ~ 
0.1 Mpc^/a^, whereas our measured A values are 100 times 
larger. 

The calculation for magnification of the LRG by fore- 
ground galaxies is analogous and gives a similar estimate. 



Here the problem is that the foreground galaxies are well 
below L* and hence are poor lenses. The volume and path 
lengths are also less favorable. 

We therefore estimate the weak lensing effects at 1-2%. 
While not important for this work, weak lensing could 
be an issue for some applications of the cross-correlation 
method, particularly for studies involving cross-correlation 
to sub-i* galaxies or to objects with favorable fc-corrections 
at high redshifts (e.g. sub- millimeter galaxies). 
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17.1 


1.0 


0.20 <z< 0.36 


1.0 




-21.2 


6.42(15) 


0.090(9) 


7.0 


1.0 


0.20 < z < 0.28 


1.0 




-21.2 


6.72(22) 


0.089(12) 


11.1 


1.0 


0.28 <z< 0.36 


1.0 




-21.2 


6.22(18) 


0.090(12) 


7.0 


1.0 


0.20 < z < 0.36 


1.0 


9-i 


-21.2 


3.17(7) 


0.099(8) 


11.0 


1.0 


0.20 < z < 0.28 


1.0 


fl - i 


-21.2 


3.39(10) 


0.106(11) 


5.0 


1.0 


n.2S < < 0.:-!(i 


1.0 


<l - i 


-21,2 


:-!.0:H8) 


n.O<):-!(ll) 


9.5 


2.0 


0.20 <z< 0.36 


1.0 




-21.2 


15.1(4) 


0.058(10) 


14.1 


2.0 


0.20 < 2 < 0.28 


1.0 




-21.2 


15.1(6) 


0.084(14) 


7.5 


2.0 


0.28 <z< 0.36 


1.0 




-21.2 


15.2(6) 


0.046(13) 


13.7 


2.0 


0.20 < 2 < 0.36 


1.0 


9-i 


-21.2 


6.64(14) 


0.059(8) 


10.4 


2.0 


0.20 < 2 < 0.28 


1.0 


9-i 


-21.2 


7.03(22) 


0.068(12) 


11.2 


2.0 


0.28 <z< 0.36 


1.0 


9-i 


-21.2 


6.45(18) 


0.053(10) 


6.5 


4.0 


0.20 < z < 0.36 


1.0 




-21.2 


36.4(15) 


0.047(13) 


8.0 


4.0 


0.20 <z < 0.28 


1.0 




-21.2 


37.4(20) 


0.082(18) 


10.1 


4.0 


0.28 < 2 < 0.36 


1.0 




-21.2 


35.7(18) 


0.024(16) 


7.7 


4.0 


0.20 <z < 0.36 


1.0 


a - i 


-21.2 


15.3(4) 


0.041(9) 


7.3 


4.0 


0.20 <z< 0.28 


1.0 


9 - i 


-21.2 


17.4(6) 


0.061(12) 


12.2 


4.0 


0.28 < 2 < 0.36 


1.0 


9-i 


-21.2 


14.0(5) 


0.022(12) 


6.6 



NOTES.—" The rcdshift range of tlie LKG sample. 

b The minimum luminosity of the imaging galaxies in magnitudes below M* . The maximum luminosity is always M* —0.5. This table presents 
the "1.0" sample. 

Marked with g — i M imaging galaxies have been required to be within 5(g — i) mag in observed g — i color of the red sequence at the LRG 
redshift. 

The minimum Mg of the LRGs used in the fit. 
" The value of the best-fit linear relation of (poVA versus (L/ L*)^'^' , evaluated at 3.3L* for the Mg < —21.2 sample and 4.5L* for the 
Mg < —21.8 sample. These values were picked so that the errors on the slope and intercept of the best-fit line are nearly independent. The 
errors on the last digits are given in parenthesis. We use M* = —20.35. 

^ The slope of the best-fit line and its error, both divided by the value in the previous column to give a reasonable normalization (but one that 
is different between this table and Table 3). 

9 of the data with respect to the best-fit line. Samples with Mg < —21.2 have 15 luminosity bins and hence 13 degrees of freedom. Samples 
with Mg < —21.8 have 9 luminosity bins and hence 7 degrees of freedom. 
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Table 3 



l-'ns lo ()\1':hi)I':ksi 1 ii:s \i'Hsrs Lriiixosm" 



a {h~^ Mpc) 


Redshift 




Color 


Llrg 




Slope 




0.125 


0.20 < z< 0.44 


0.4 




-21.8 


0.637(13) 


0.065(7) 


7.9 


0.125 


0.20 <z< 0.28 


0.4 




-21.8 


0.645(27) 


0.085(18) 


2.7 


0.125 


0.28 <z< 0.36 


0.4 




-21.8 


0.673(20) 


0.063(12) 


7.2 


0.125 


0.36 <z< 0.44 


0.4 




-21.8 


0.616(20) 


0.068(11) 


5.1 


0.125 


0.20 <z< 0.44 


0.4 


9-i 


-21.8 


0.468(8) 


0.071(7) 


4.7 


0.125 


0.20 < z < 0.28 


0.4 


9-i 


—21.8 


0.530(21) 


0.093(17) 


3.3 


0.125 


0.28 < z < 0.36 


0.4 




—21.8 


0.481(13) 


0.072(11) 


2.9 


0.125 


0.36 < z < 0.44 


0.4 




—21.8 


0.433(11) 


0.068(10) 


2.5 


0.25 


0.20 < z < 0.44 


0.4 




-21.8 


1.314(23) 


0.084(7) 


6.0 


0.25 


0.20 <z < 0.28 


0.4 




-21.8 


1.45(5) 


0.091(15) 


3.8 


0.25 


0.28 < 2 < 0.36 


0.4 




-21.8 


1.355(38) 


0.067(11) 


3.5 


0.25 


0.36 < z < 0.44 


0.4 




-21.8 


1.219(35) 


0.098(11) 


3.3 


0.25 


0.20 < z < 0.44 


0.4 




-21.8 


0.870(15) 


0.090(7) 


7.9 


0.25 


0.20 < z < 0.28 


0.4 


9-i 


—21.8 


1.013(35) 


0.103(15) 


3.4 


0.25 


0.28 < z < 0.36 


0.4 


n — i 


—21.8 


0.911(24) 


0.088(10) 


5.6 


0.25 


36 ^ 2 ^ 44 


0.4 


9 ~ 


—21.8 


782^21^ 


084no^ 


7.2 


0.5 


0.20 <z< 0.44 


0.4 




-21.8 


2.21(5) 


0.083(8) 


4.8 


0.5 


0.20 < 2 < 0.28 


0.4 




-21.8 


2.47(9) 


0.092(17) 


0.8 


0.5 


0.28 < 2 < 0.36 


0.4 




-21.8 


2.37(8) 


0.075(14) 


5.2 


0.5 


0.36 < 2 < 0.44 


0.4 




-21.8 


2.00(8) 


0.091(14) 


5.3 


0.5 


0.20 < 2 < 0.44 


0.4 


9 - i 


-21.8 


1.342(28) 


0.098(8) 


10.4 


0.5 


0.20 < 2 < 0.28 


0.4 


9-i 


—21.8 


1.63(6) 


0.106(17) 


3.1 


0.5 


0.28 < 2 < 0.36 


0.4 


(7 — i 
y " 


-21.8 


1 41(5) 


0.103(13) 


10.0 


0.5 


0.36 < 2 < 0.44 


0.4 


o — i 
y " 


-21.8 


1 18('4~) 


0.097(11) 


10.7 


1.0 


0.20 < 2 < 0.44 


0.4 




-21.8 


3.68(12) 


0.067(11) 


12.2 


1.0 


0.20 < 2 < 0.28 


0.4 




-21.8 


4.00(25) 


0.075(21) 


4.9 


1.0 


0.28 < 2 < 0.36 


0.4 




-21.8 


3.92(19) 


0.080(20) 


4.6 


1.0 


0.36 <z< 0.44 


0.4 




-21.8 


3.39(18) 


0.043(18) 


12.1 


1.0 


0.20 < z < 0.44 


0.4 


9-i 


-21.8 


1.98(5) 


0.068(9) 


8.2 


1.0 


0.20 < z < 0.28 


0.4 


y *' 


—21.8 


2.34(11) 


0.076(19) 


1.8 


1.0 


0.28 < 2 < 0.36 


0.4 


9-i 


-21.8 


2.09(9) 


0.098(19) 


3.4 


1.0 


0.36 < 2 < 0.44 


0.4 


9-i 


-21.8 


1.78(7) 


0.044(12) 


4.2 


2.0 


0.20 < 2 < 0.44 


0.4 




-21.8 


8.19(33) 


0.064(12) 


11.9 


2.0 


0.20 < 2 < 0.28 


0.4 




-21.8 


8.7(6) 


0.078(24) 


6.9 


2.0 


0.28 < 2 < 0.36 


0.4 




-21.8 


7.9(5) 


0.076(26) 


9.2 


2.0 


0.36 < 2 < 0.44 


0.4 




-21.8 


8.1(5) 


0.056(20) 


8.2 


2.0 


0.20 < 2 < 0.44 


0.4 


9 - i 


-21.8 


4.07(11) 


0.055(9) 


10.1 


2.0 


0.20 < 2 < 0.28 


0.4 


9-i 


-21.8 


4.57(24) 


0.069(21) 


5.8 


2.0 


0.28 < 2 < 0.36 


0.4 


9-i 


-21.8 


4.19(19) 


0.062(18) 


4.2 


2.0 


0.36 <z< 0.44 


0.4 


9-i 


-21.8 


3.77(16) 


0.040(14) 


4.2 


4.0 


0.20 < z < 0.44 


0.4 




-21.8 


20.2(13) 


0.031(17) 


5.2 


4.0 


0.20 <z< 0.28 


0.4 




-21.8 


24.3(19) 


0.061(27) 


3.5 


4.0 


0.28 <z< 0.36 


0.4 




-21.8 


20.5(20) 


-0.003(28) 


5.0 


4.0 


0.36 < 2 < 0.44 


0.4 




-21.8 


18.4(17) 


0.035(27) 


6.6 


4.0 


0.20 < 2 < 0.44 


0.4 


9-i 


-21.8 


9.29(35) 


0.027(11) 


5.4 


4.0 


0.20 < 2 < 0.28 


0.4 


9-i 


-21.8 


12.0(7) 


0.055(20) 


3.4 


4.0 


0.28 < 2 < 0.36 


0.4 


9-i 


-21.8 


9.0(5) 


0.009(19) 


3.8 


4.0 


0.36 < 2 < 0.44 


0.4 


9-i 


-21.8 


8.5(5) 


0.021(18) 


5.8 



NOTES. — As Table 2, but for the higher luminosity "0.4" imaging sample, which allows us to use higher redshift LRGs if we restrict to high 
luminosity LRGs. 



